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ABSTRACT 


During  the  final  year  of  our  project,  we  have  calibrated  the  use  of  surface-wave  magnitudes  (M^),  measured  at 
regional  distances  as  a  rapid  and  robust  estimator  of  seismic  moment.  We  have  used  the  Russell  (2006)  variable- 
period  surface-wave  magnitude  formula  to  convert  to  seismic-moment  magnitude,  M^,  at  local  to  regional 
distances  using  global  datasets.  In  this  pilot  study,  the  Russell  technology  was  applied  to  169  North  American 
events  with  3.2  <  <  6.5  at  distances  ranging  from  48  to  5,268  km.  The  technique  uses  a  time-domain  magnitude 

estimation  procedure  that  employs  zero-phase  Butterworth  filters  to  effectively  measure  Rayleigh-wave  Airy  phase 
amplitudes  at  local  and  near-regional  distances  (e.g.  ,<  1000  km).  This  allows  for  surface  wave  magnitudes  to  be 
estimated  within  minutes  of  the  initiation  of  a  seismic  event.  Of  the  7,370  event-station  pairs,  more  than  half  (4,051) 
of  the  measurements  were  at  distances  <  1 ,000  km. 

The  Ms  estimates  were  regressed  against  moment  magnitudes  (M^)  estimated  from  7^-wave  modeling  and/or 
Rayleigh-  and  Love-wave  spectral  amplitudes  (Herrmann  et  al.,  2008).  can  be  estimated  using  the  relationship: 

=  1.951  +  0.649  Ms.  The  observed  scatter  in  the  estimated  was  approximately  ±0.2  magnitude  units.  The 
residuals  between  true  and  M^-predicted  M^  have  a  definable  faulting  mechanism  effect,  especially  when  strike-slip 
events  are  compared  to  those  with  other  mechanisms.  Preliminary  results  suggest  that  our  Ms:M^  relationship  for 
North  America  is  also  transportable  to  the  Middle  East  and  the  Korean  Peninsula. 

We  have  also  determined  a  methodology  to  estimate  the  Ms  detection  thresholds  for  the  Russell  formula.  Broadband 
noise  estimates  for  any  seismic  station,  which  are  typically  in  units  of  acceleration  m^/s^  and  decibels,  can  be 
converted  to  nanometers  (nm)  and  input  into  the  Russell  equation  for  variable-period  surface  waves.  We  propagate 
these  noise  estimates  at  periods  (T)  between  8  and  40  s  to  distances  (A)  between  1  and  40  degrees.  The  Ms{T,^  were 
estimated  and  contoured  providing  the  optimal  detection  threshold  for  each  station.  We  have  estimated  the  Ms(T,S) 
for  the  stations  of  the  Global  Seismographic  Network  (GSN)  using  broadband  noise  estimates  supplied  by  Berger  et 
al.  (2004). 
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OBJECTIVES 

Developing  a  methodology  for  calculating  surface  wave  magnitudes  that  is  valid  at  both  regional  and  teleseismic 
distances,  applicable  to  events  of  variable  sizes  and  signal-to-noise  ratios,  calibrated  for  variable  structure  and 
propagation,  and  easy  to  automate  in  an  operational  setting,  is  an  important  monitoring  goal.  Our  objectives  are  to 
create  such  a  methodology,  and  to  use  it  to  lower  Ms  estimation  and  detection  thresholds.  We  hope  that  the  method 
will  provide  a  seamless  tie  between  Ms  estimation  at  regional  and  teleseismic  distances. 

To  accomplish  our  objectives,  we  developed  the  Matlab  program  EVALSURF  (Bonner  et  al.,  2006),  which 
estimates  variable-period  (8  <  T  <  25  s;  recently  updated  to  40  s)  Rayleigh-wave  magnitudes  using  the  Russell 
(2006)  and  Ms(VMAX)  measurement  technique  (Bonner  et  al.,  2006)  for  comparison  to  the  historical  formulas 
of  Marshall  and  Basham  (1972)  and  Rezapour  and  Pearce  (1998).  The  program  uses  the  updated  Lawrence 
Livermore  National  Laboratory  group  velocity  models  for  Eurasia  (Pasyanos,  2005)  to  identify,  phase  match 
filter,  and  extract  the  fundamental-mode  Rayleigh  waves  for  analysis.  During  the  past  year,  we  calibrated  the  use 
of  surface  wave  magnitudes  (M^),  measured  at  regional  distances  as  a  rapid  and  robust  estimator  of  seismic 
moment.  We  have  also  determined  a  methodology  to  estimate  the  Ms  detection  thresholds  for  European  and 
Asian  Global  Seismographic  Network  (GSN)  stations. 

RESEARCH  ACCOMPLISHED 
Surface  Wave  Magnitudes 

Many  surface  wave  magnitude  scales  have  been  developed  during  the  past  century.  They  have  either  been  based  on 
empirical  (Gutenberg,  1945;  Vanek  et  al.,  1962;  Marshall  and  Basham,  1972)  or  theoretical  (Rezapour  and  Pearce, 
1998;  Stevens  and  McLaughlin,  2001)  aspects  of  Rayleigh- wave  propagation.  Knowing  how  and  when  to  apply 
these  formulas  has  often  been  confusing,  as  some  of  these  formulas  were  developed  strictly  for  ~20-second  surface 
waves  at  teleseismic  distances.  Others  were  developed  for  regional  and  teleseismic  distances  and  variable  periods 
(>  10s).  What  was  needed  was  a  seamless  relationship  between  estimating  Ms  at  regional  and  teleseismic  distances 
for  events  of  a  wider  range  of  magnitudes. 

Russell  (2006)  addressed  this  need  by  developing  a  time-domain  method  for  measuring  surface  waves  with  minimal 
digital  processing,  using  zero-phase  Butterworth  filters.  The  method  can  effectively  measure  surface-wave 
magnitudes  at  both  regional  and  teleseismic  distances,  at  variable  periods  between  8  and  25  s.  The  magnitude 
equation  is 


M 


s(b) 


=  logfaj  +  ilog(sin(zl))+0.003/^ 


I  T 


zl-0.661og|^ 


-log(/J-0.43, 


(1) 


where  a\^  is  the  amplitude  of  the  Butterworth-filtered  surface  waves  (zero-to-peak  in  nanometers)  and  fc  is  the  filter 
frequency  of  a  zero-phase  Butterworth  band-pass  filter  with  comer  frequencies  1/T -f^  1/T +fc.  At  the  reference 
period  r=20  s,  the  equation  is  equivalent  to  Von  Seggem's  formula  (1977)  scaled  to  Vanek  et  al.  (1962)  at  50 
degrees.  For  periods  8<r<25,  the  equation  is  corrected  to  r=20  sec,  accounting  for  source  effects,  attenuation,  and 
dispersion.  This  technique  is  often  referred  to  as  Ms(VMAX)  for  Variable-period,  MAXimum  amplitude  magnitude 
estimates  (Bonner  et  al.,  2006). 

The  advantages  of  the  M^VMAX)  technique  include 

•  The  technique  allows  for  time-domain  measurements  of  surface-wave  amplitudes,  giving  an  analyst  the 
ability  to  easily  and  visually  confirm  that  the  pick  is  correct  and  is  an  actual  Rayleigh  wave. 

•  It  allows  for  surface-wave  magnitudes  to  be  measured  at  local  and  regional  distances  where  traditional 
20-s  magnitudes  cannot  be  used.  These  magnitudes  are  not  biased  with  respect  to  teleseismic  estimates. 

•  The  application  of  narrow-band  Butterworth-filtering  techniques  appropriately  handles  the  Airy  phase. 

•  The  method  is  variable  period  and  not  restricted  to  near  20-s  period,  and  the  analyst  is  allowed  to  measure 
the  surface  wave  magnitude  where  the  signal  is  largest. 
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M,(VMAX)  Methodology 

To  estimate  ^/^(VMAX),  we  apply  a  series  of  zero-phase  3^^-order  Butterworth  filters  to  the  data  with  the  comer 
frequencies  (l/T)-fc,  {l/T)+fc^  where 


^  Gmin 

“  t4a  ' 


(2) 


For  continental  paths  between  8  and  25  s,  Gmin,  which  is  based  on  dispersion  effects,  is  set  to  0.6.  The  center 
periods  are  placed  at  1 -second  intervals  between  8  and  25  s.  We  constmct  the  envelope  function  of  the  filtered  signal 
and  measure  the  maximum  zero-to-peak  amplitude  in  a  group  velocity  window  between  2.0  and  4.0  km/s.  An 
analyst  then  visually  confirms  that  the  correct  waveform  feature  is  being  measured — a  benefit  of  using  a 
time-domain  measurement.  The  procedure  can  be  automated  based  on  signal-to-noise  requirements  within  the 
window.  In  Figure  1  we  show  examples  of  filter  panels  from  an  Mvi;=4.83  Montana  earthquake  recorded  at  station 
DUG  at  a  near-regional  distance  of  509  km.  Note  the  differences  in  the  Rayleigh-wave  amplitudes  between  20-s 
periods,  where  traditional  My  is  measured,  and  8-14  s  periods. 

We  record  the  maximum  amplitude  in  each  of  the  18  filter  bands  and  then  use  Equation  1  to  calculate  a  variable- 
period  surface-wave  magnitude.  As  a  result,  18  different  magnitudes  are  estimated  for  each  station  recording  the 
event  (Figure  2).  We  then  search  the  variable-period  filtered  data  (Figure  1)  to  determine  the  period  of  the  maximum 
Airy-corrected  amplitude.  We  use  the  magnitude  calculated  at  that  period  as  the  final  estimate.  As  an  example,  the 
period  of  maximum  amplitude  for  the  DUG  recording  occurred  at  T=10  s.  We  note  that  for  this  Montana  event,  all 
of  the  M^VMAX)  estimates  were  made  at  periods  of  less  than  18  sec.  We  do  the  same  processing  for  all  stations 
that  recorded  the  event  and  then  estimate  a  network  M^(VMAX)  (4.42),  the  interstation  standard  deviation  (0.13), 
and  surface-wave  magnitude  “noise”  floors  at  each  station  (see  Figure  2). 
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Figure  1.  Filter  combs  of  Rayleigh  waves  at  periods  of  8,  10,  12, ..,24  seconds  for  a  Montana  earthquake 
recorded  at  DUG.  The  first  blue  line  is  4.0  km/s,  while  the  second  is  2.0  km/s.  The  red  line  marks 
the  largest  amplitude.  Amplitudes  are  in  nm. 
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Figure  2.  M5(VMAX)  results  for  the  Rayleigh  waves  recorded  from  a  Montana  earthquake.  The  solid  lines 
represent  signals  while  the  dashed  lines  are  pre-surface-wave  noise  levels.  The  station  magnitudes 
were  averaged  to  form  a  network  Ms(VMAX)  of  4.42  for  this  M^=4.83  event. 


Application  in  North  America 

We  have  used  the  M^VMAX)  method  to  estimate  surface-wave  magnitudes  for  169  events  in  North  America.  The 
events  were  part  of  a  database  of  North  American  events  for  which  Herrmann  et  al.  (2008)  measured  anM^,  using 
either  broadband  waveform  modeling  and/or  Rayleigh-  and  Love-wave  spectral  amplitudes.  The  events  include  most 
types  of  faulting  focal  mechanisms  (Zoback,  1992).  Figure  3  shows  the  geographic  distribution  of  these  events, 
which  were  3.2  <  <6.5,  and  stations.  My(VMAX)  was  estimated  for  stations  at  distances  ranging  from  48  to 

5,268  km.  Of  the  7,370  event-station  pairs,  4,051  of  the  measurements  were  at  distances  <  1,000  km. 


Figure  3.  Earthquakes  (yellow  stars)  and  stations  (red  circles)  used  in  a  pilot  study  to  examine  the 
effectiveness  of  using  regional  surface-wave  magnitudes  to  estimate  M^. 
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An  example  of  the  processing  for  the  21  February  2008  Wells,  Nevada  (M^=5.88)  event  is  shown  in  Figure  4.  This 
event  was  located  within  the  heart  of  Earthscope’s  Transportable  Array,  and  thus  was  recorded  on  hundreds  of 
stations  within  North  America  alone.  Because  the  Ms(YMAX)  technology  can  estimate  unbiased  surface  wave 
magnitudes  at  epicentral  distances  as  close  as  50  km,  we  generated  506  surface-wave  magnitude  estimates  for  this 
event  with  only  a  very  small  azimuthal  gap.  Even  though  the  focal  mechanism  for  this  normal  fault  event  shows 
differences  in  the  surface  wave  radiation  patterns  at  10,  15,  20,  and  25  s  (Herrmann  et  al.,  2008),  the  radiation 
pattern  of  the  estimates  shows  no  significant  radiation  patterns.  This  is  due  in  part  to  the  fact  that  the  variable- 
period  methodology  picks  the  surface-wave  estimate  at  the  largest  amplitude  and  often  migrates  away  from  periods 
with  radiation  pattern  nulls.  We  determined  an  Ms(VMAX)  of  5.99  for  the  event. 
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Figure  4.  M^VMAX)  estimation  for  the  21  February  2008  Wells,  NV  earthquake.  Shown  are  the  stations 
used  to  estimate  the  variable-period  surface-wave  magnitude,  the  magnitude  analysis  at  periods 
between  8-25  seconds,  surface  wave  radiation  patterns  from  Herrmann  et  al.  (2008),  and  the  Mg 
estimates  as  a  function  of  event-to-station  azimuth. 

Relationship  Between  and  M^VMAX) 

We  derived  a  relationship  between  and  M^VMAX)  using  the  169  North  American  earthquakes.  The  Mg 
estimates  were  regressed  (Figure  5)  against  moment  magnitudes  (M^)  and  log  estimated  from  E’-wave  modeling 
and/or  Rayleigh-  and  Love-wave  spectral  amplitudes.  can  be  estimated  using  the  relationship:  =  1.951  + 

0.649  Mg.  To  demonstrate  the  observed  scatter  in  the  method,  we  show  the  residual  Mg  estimates  as  a  function  of 
period  and  distance  in  Figure  6.  With  the  exception  of  T=8  sec,  where  highly  variable  crustal  structure  is 
problematic,  the  residual  surface-wave  magnitude  estimates  are  almost  always  within  ±0.2  magnitude  units  of  the 
mean  value.  And  there  are  no  significant  distance  trends  observable  in  the  data  at  periods  greater  than  10  sec. 
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Figure  5.  Linear  regressions  of  (left)  and  log  Mg  (right)  versus  M5(VMAX).  The  moments  were 
determined  from  body  wave  modeling  and/or  Love  and  Rayleigh-wave  spectral  amplitudes 
(Herrmann  et  al.,  2008). 
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Figure  6.  Scatter  in  the  estimated  /^^(VMAX)  estimates  as  a  function  of  the  evaluation  period  and  distance. 

The  scatter  (Figure  7)  between  the  observed  and  predicted  shows  mechanism  dependence  for  normal  and  strike- 
slip  faulting  sources.  Each  focal  mechanism  was  classified  as  either  normal  (NF)  strike-slip  (SS),  thrust  (TF),  or 
oblique-slip  variations  of  normal  and  thrust  faults  using  the  Zoback  (1992)  classification  scheme.  For  normal  faults, 
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the  predicted  is  -0.04  magnitude  units  larger  than  the  observed  A  similar  bias  might  be  noted  for  thrust 
faulting  as  more  estimates  are  obtained  in  the  future.  On  the  contrary,  for  strike-slip  faults,  the  predicted  estimate  is 
-0.07  magnitude  units  smaller  than  the  observed  M^i;.  The  reason  for  these  differences  results  from  ourM^  estimates 
only  being  based  on  the  Rayleigh  waves.  Strike-slip  events  contain  a  significant  amount  of  their  moment  from 
energy  in  the  Love  waves,  which  are  not  considered  in  our  technique.  The  results  show  there  is  little,  if  any,  depth 
dependence  on  the  predicted  for  events  with  crustal  depths. 
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Figure  7.  The  dependence  of  the  M^-predicted  as  a  function  of  mechanism  and  depth.  The  mechanisms 

are  based  on  the  Zoback  (1992)  classiBcation  scheme.  The  residuals  between  the  observed  and 
M^-predicted  have  a  definable  faulting  mechanism  effect,  especially  when  strike-slip  events  are 
compared  to  those  with  other  mechanisms. 


Predicting  M^VMAX)  Detection  Thresholds 

Broadband  noise  estimates  can  be  used  to  estimate  the  M^YMAX)  threshold  floor  for  any  seismic  station.  Berger  et 
al.  (2004)  compiled  broadband  ambient  noise  estimates  for  the  GSN.  Peter  Davis  (Executive  Director,  Project  IDA 
at  the  University  of  California,  San  Diego)  provided  us  with  the  noise  estimates  for  the  GSN  stations.  As  an 
example,  the  left  subplot  of  Figure  8  shows  the  25^^  percentile  noise  estimate  for  station  AAK  compared  to  the  low 
noise  model  for  the  GSN. 

We  converted  these  noise  estimates,  which  were  in  units  of  acceleration  m^/s^  and  decibels,  to  nanometers  (nm)  and 
input  them  into  Equation  1  for  variable-period  surface  waves.  We  propagated  these  noise  estimates  at  periods  (T) 
between  8  and  25  s  to  distances  (A.between  1  and  60  degrees.  The  Ms(T,A)  were  estimated  and  contoured 
(right  subplot  of  Figure  8).  This  colormap  provides  the  optimal  detection  threshold  for  each  station  for  a  signal-to- 
noise  ratio  of  2.  If  the  event  has  an  Mg  smaller  than  these  values  at  a  particular  distance  or  period,  we  will  probably 
discard  it  due  to  the  inability  to  verify  that  it  is  a  Rayleigh  wave  because  of  the  "noise-based"  Ms  background  level. 
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Using  the  GSN  noise  estimates  (examples  are  shown  in  Figures  8  and  9),  we  determine  the  maximum  distance  in 
which  we  could  record  an  M^=2.5  and  My=3.0.  This  distance  is  plotted  as  the  radius  of  a  circle  centered  on  the 
station  in  Figures  10  and  1 1  for  Europe  and  Asia.  The  blue  circle  shows  the  detection  limits  for  the  M^VMAX) 
technique,  which  considers  periods  between  8  and  25  s.  The  red  circle  shows  detection  distances  if  only  periods 
between  17  and  23  s  were  considered,  which  is  the  traditional  bandwidth  for  surface  wave  magnitude  estimation.  For 
Ms=2.5,  the  M^YMAX)  technique  has  a  mean  detection  distance  of  9°  compared  with  7°  for  traditional 
17-23-second  measurements.  For  Ms=3.0,  the  mean  detection  distance  is  26°  for  M^YMAX)  and  22°  for  the 
traditional  bandwidth. 
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Figure  8.  M5(VMAX)  detection  threshold  for  station  ABKT.  Left:  Long-period  noise  estimates  (blue  line)  for 
the  station  compared  to  the  GSN  low  noise  model  (solid  black  line).  Right:  Detection  magnitude 
thresholds  contoured  as  a  function  of  distance  and  period. 
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Figure  9.  M^(VMAX)  detection  threshold  for  station  RAYN.  Left:  Long-period  noise  estimates  (blue  line)  for 
the  station  compared  to  the  GSN  low-noise  model  (solid  black  line).  Right:  Detection  magnitude 
thresholds  contoured  as  a  function  of  distance  and  period. 
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Figure  10.  Detection  thresholds  for  Europe  and  Asian  GSN  stations  for  an  M^(VMAX)=2.5.  The  blue  curve 
shows  the  detection  threshold  of  Ms(VMAX)  between  8-  and  25-second  periods.  The  red  curve 
shows  the  detection  thresholds  only  using  17-22-second  surface  waves,  which  is  the  traditional 
bandwidth  for  analysis. 
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Figure  11.  Detection  thresholds  for  Europe  and  Asia  for  an  M5(VMAX)=3.  0. 
CONCLUSIONS  AND  RECOMMENDATIONS 


The  variable-period  surface  wave  magnitude  M^VMAX)  continues  to  provide  a  better  and  more  consistent  estimate 
of  source  size,  particularly  for  smaller  events  and  at  shorter  distances.  During  the  past  year,  we  applied  Ms(VMAX) 
to  169  North  American  events  with  3.2  <  <  6.5  at  distances  ranging  from  48  to  5,268  km.  The  Ms  estimates  were 

regressed  against  moment  magnitudes  (My^)  estimated  from  T’-wave  modeling  and/or  Rayleigh-  and  Love-wave 
spectral  amplitudes.  The  resulting  Ms.'My^  relationship  has  scatter  on  the  order  of  ±0.2  magnitude  units  in  the 
estimated  My^.  The  residuals  between  true  and  My-predicted  Myy  have  a  definable  faulting  mechanism  effect, 
especially  when  strike-slip  events  are  compared  to  those  with  other  mechanisms.  This  could  have  important 
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ramifications  for  the  Ms.-nih  discriminant  in  regions  where  a  particular  mechanism  dominates  the  seismicity.  This 
exercise  provides  a  tool  for  rapid  estimation  of  events  with  >3.5  for  North  America.  Preliminary  application 
of  the  same  procedure  to  two  other  regions  suggests  that  this  Ms:M^  relationship  is  also  transportable  to  the  Middle 
East  and  Korean  Peninsula  (Figure  12). 


Figure  12.  Comparison  of  log  Mo  versus 
M5(VMAX)  for  earthquakes  in 
North  America  (green),  Korea  (red) 
and  the  Middle  East  (black).  Note 
that  the  Korean  earthquakes  are 
mostly  strike-slip. 
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